Search for: All records

Our knowledge of stellar evolution is driven by one-dimensional (1D) simulations. 1D models, however, are severely limited by uncertainties on the exact behaviour of many multidimensional phenomena occurring inside stars, affecting their structure and evolution. Recent advances in computing resources have allowed small sections of a star to be reproduced with multi-D hydrodynamic models, with an unprecedented degree of detail and realism. In this work, we present a set of 3D simulations of a convective neon-burning shell in a 20 M⊙ star run for the first time continuously from its early development through to complete fuel exhaustion, using unaltered input conditions from a 321D-guided 1D stellar model. These simulations help answer some open questions in stellar physics. In particular, they show that convective regions do not grow indefinitely due to entrainment of fresh material, but fuel consumption prevails over entrainment, so when fuel is exhausted convection also starts decaying. Our results show convergence between the multi-D simulations and the new 321D-guided 1D model, concerning the amount of convective boundary mixing to include in stellar models. The size of the convective zones in a star strongly affects its structure and evolution; thus, revising their modelling in 1D will have important implications for the life and fate of stars. This will thus affect theoretical predictions related to nucleosynthesis, supernova explosions, and compact remnants.

 

We present the state-of-the-art single-zone nuclear reaction networkWinNet, which is capable of calculating the nucleosynthetic yields of a large variety of astrophysical environments and conditions. This ranges from the calculation of the primordial nucleosynthesis, where only a few nuclei are considered, to the ejecta of neutron star mergers with several thousands of involved nuclei. Here we describe the underlying physics and implementation details of the reaction network. We additionally present the numerical implementation of two different integration methods, the implicit Euler method and Gears method, along with their advantages and disadvantages. We furthermore describe basic example cases of thermodynamic conditions that we provide together with the network and demonstrate the reliability of the code by using simple test cases. With this publication,WinNetwill be publicly available and open source at GitHub and Zenodo.

 

ABSTRACT The origin of the proton-rich trans-iron isotopes in the Solar system is still uncertain. Single-degenerate thermonuclear supernovae (SNIa) with n-capture nucleosynthesis seeds assembled in the external layers of the progenitor’s rapidly accreting white dwarf (RAWD) phase may produce these isotopes. We calculate the stellar structure of the accretion phase of five white dwarf (WD) models with initial masses ≥ 0.85 $\, \mathrm{M}_\odot$ using the stellar code mesa The near-surface layers of the 1, 1.26, 1.32 and 1.38 $\, \mathrm{M}_\odot$ models are most representative of the regions in which the bulk of the p nuclei are produced during SNIa explosions, and for these models we also calculate the neutron-capture nucleosynthesis in the external layers. Contrary to previous RAWD models at lower mass, we find that the H-shell flashes are the main site of n-capture nucleosynthesis. We find high neutron densities up to several 1015 cm−3 in the most massive WDs. Through the recurrence of the H-shell flashes these intermediate neutron densities can be sustained effectively for a long time leading to high-neutron exposures with a strong production up to Pb. Both the neutron density and the neutron exposure increase with increasing the mass of the accreting WD. Finally, the SNIa nucleosynthesis is calculated using the obtained abundances as seeds. We obtain solar to supersolar abundances for p-nuclei with A > 96. Our models show that SNIa are a viable p-process production site. 

Meteoritic analysis demonstrates that radioactive nuclei heavier than iron were present in the early Solar System. Among them, 129I and 247Cm both have a rapid neutron-capture process (r process) origin and decay on the same timescale (≃ 15.6 Myr). We show that the 129I/247Cm abundance ratio in the early Solar System (438±184) is immune to galactic evolution uncertainties and represents the first direct observational constraint for the properties of the last r-process event that polluted the pre-solar nebula. We investigate the physical conditions of this event using nucleosynthesis calculations and demonstrate that moderately neutron-rich ejecta can produce the observed ratio. We conclude that a dominant contribution by exceedingly neutron-rich ejecta is highly disfavoured.